Strain controlled oxygen vacancy formation and ordering in CaMn0 3 
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We use first-principles calculations to investigate the stability of bi-axially strained Pnma per- 
ovskite CaMnOs towards the formation of oxygen vacancies. Our motivation is provided by promis- 
ing indications that novel material properties can be engineered by application of strain through 
coherent heteroepitaxy in thin films. While it is usually assumed that such epitaxial strain is accom- 
modated primarily by changes in intrinsic lattice constants, point defect formation is also a likely 
strain relaxation mechanism. This is particularly true at the large strain magnitudes (>4%) which 
first-principles calculations often suggest are required to induce new functionalities. We find a strong 
dependence of oxygen vacancy defect formation energy on strain, with tensile strain lowering the 
formation energy consistent with the increasing molar volume with increasing oxygen deficiency. In 
addition, we find that strain differentiates the formation energy for different lattice sites, suggesting 
its use as a route to engineering vacancy ordering in epitaxial thin films. 



I. INTRODUCTION 

Materials of the ABO3 perovskite family and derived 
structures exhibit a large variety of functional proper- 
ties which are of fundamental and technological inter- 
est. Intrinsic and extrinsic point defects, such as oxy- 
gen vacancies or cation non-stoichiometry always ex- 
ist in these materials and are in some cases desirable 
or even crucial for their properties. For example oxy- 
gen vacancies enable ionic conductivity in high tem- 
perature electrochemial devices such as oxygen sensors, 
SOFC cathodes and mixed conducting membranes^'. 
Deviations from ideal stoichiometry and charge balance 
are essential for introducing the carriers required for 
conventional- or super-conductivity, or carrier-mediated 
ferromagnetic ordering. The arrangement of defects can 
also be important, with reports of ordered oxygen vacan- 
cies inducing high-temperature superconductivity^ and 
multiferroicitySl. The electronic structure of interfaces 
in oxide heterostructures such as the formation of two- 
dimensional electron gases (2DEG)0, possibly relies on 
cation interdiffusion and will thus be critically affected by 
the concentration of cation vacancies required for cation 
migration. Point defects can also be detrimental, how- 
ever, particularly if insulating behavior is desired, such 
as in ferroelectric or capacitor applications, or in achiev- 
ing intrinsic antiferromagnetic ordering through superex- 
change. 

In bulk perovskite oxides the presence of oxygen va- 
cancies is well established to result in an expansion of 
the crystal lattice due to underbonding caused by the 
two extra electrons from the missing oxygen located in 
non-bonding orbitals^l This phenomenon, known as 
chemical expansion, is most pronounced for perovskites 
containing transition metal ions that can adopt a vari- 
ety of valence states and hence can readily accommodate 
the associated change in formal charge^. Prominent ex- 
amples are acceptor doped LaMn03 and donor doped 
SrTiOi^LUJ where significant volume expansions are ob- 



served function of the oxygen partial pressure. 

It has been demonstrated during the last decade that 
functional behavior in perovskites can be enabled or en- 
hanced when they are prepared as thin films. For ex- 
ample, ferroelectricity has been reported in thin films 
of SrTi03, which is a quantum paraelectric in its bulk 
form, and enhanced polarization has been observed in 
thin films of ferroelectric BaTiOgpESl Reports of cou- 
pling magnetic and ferroelectric order p aram eters to yield 
thin-film magnetoelectric multiferroic d 14 l 1 5 I have spiked 
tremendous interest. The improved behavior is usu- 
ally attributed to the bi-axial epitaxial strain introduced 
when the film is grown coherently on a substrate with 
a different lattice constant. In turn, it is assumed that 
the change in in-plane lattice constant associated with 
epitaxial strain is accommodated through structural dis- 
tortions such as changing the internal bond lengths (top 
row of Fig. [T]), or changing the angle or pattern of ro- 
tations and tilts of the oxygen octahedra (center row of 
Fig. [TpSHm Indeed first-principles electronic-structure 
calculations within this approximation have proven suc- 
cessful in predi cting an d explaining many of the observed 
novel behavior d 14 * 19 * 20 ^ for a review see RefAS. 

In addition to changes in internal bond lengths or octa- 
hedral rotations, another possible effect of the change in 
lattice parameter forced by coherent epitaxy is a change 
in the defect profile. Such a response is likely because 
both strain and point defects affect the molar volume 
and hence the Gibbs free energy resulting in strain de- 
pendent defect formation energies. The increase in vol- 
ume associated with biaxial tensile strain is likely to in- 
crease the concentration of oxygen vacancies (Fig. [l] 
lower right panel). In the opposite direction, compressive 
strain should increase the concentration of cation vacan- 
cies (Fig. [l] lower left panel) due to the concomitant 
charge balancing oxidation of cations with smaller ionic 
radii at higher oxidation states. Coherent epitaxial strain 
could thus be used as a parameter to control point defect 
populations. Conversely, the spontaneous formation of 
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FIG. 1. Schematic illustration of possible mechanisms for epi- 
taxial strain accommodation in perovskites. Only B-site octa- 
hedra are shown, A-site cations being omitted for clarity. Top 
row: accommodation by bond-length changes, center row: ac- 
commodation by changes in octahedral tilt (a) and rotation 
(/?) magnitudes, bottom row: accommodation by point defect 
formation. In the bottom right we illustrate the chemical ex- 
pansion associated with oxygen vacancy formation, which is 
the subject of this work; the bottom left shows the "chemical 
contraction" caused by cation vacancy formation with con- 
comitant oxidation of cations^. 

oxygen vacancies could act as a route of chemistry me- 
diated strain relaxatio n 1 10 * 21 ! and allow large lattice mis- 
matches to be accommodated in heterostructures. We 
explore these possibilities here. 

We choose the perovskite CaMnC>3 as our model com- 
pound for analyzing the strain dependence of the oxy- 
gen vacancy (Vo) formation reaction. CaMn03 forms in 
the widely adopted Pnma modification of the ideal per- 
ovskite structure (Fig. |2j, and this tilt pattern has been 
shown to remain stable throughout t he entir e range of 
likely experimentally accessible 

strain JT5l20l22l 

The oxy- 
gen vacancy has a charge of +2 relative to the perfect 
lattice, usually written as V" in Kroger- Vink notation, 
but here Vo is used for simplicity. This relative charge of 
Vo has to be compensated by point defects with oppo- 
site charge to retain charge neutrality. The compensating 
point defect is Mn 3+ on Mn 4+ sites with a relative charge 
of -1, Mn Mn in Kroger- Vink notation. The chemical re- 
action that we study can be written as 

2Mn 4+ + 2 - ^ 2Mn 3+ + V + V20 2 (g) (1) 

and has an experimental enthalpy of 1.89±0.04 eV in 
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FIG. 2. The Pnma crystal structure of bulk, unstrained 
CaMnOs. The two inequivalent Vo positions, which are ad- 
jacent to Mn ions in the strained ac-plane (IP) and perpen- 
dicular to it (OP) are indicated. Color code: Ca: blue, Mn: 
purple, O: red. 

bulk CaMnOjpS] Since high-spin Mn 3 + (0.645 A) is 
larger than Mn 4+ (0.53 Aj22 it is clear that reaction 
([!} should be accompanied by a considerable volume ex- 
pansion, which indeed is supported by experiment; J 25 ' 26 4 
Since volume increases (decreases) with tensile (compres- 
sive) strain, we expect the Vo formation energy to be sen- 
sitive to strain. In addition to this volume effect, strain 
may also affect the Vo formation energy through its ef- 
fect on the electronic energy levels. The electrons from 
the loss of O2 which reduce Mn 4+ to Mn 3+ occupy the 
e g orbitals, the energies of which are strongly affected by 
the coordination geometry of the surrounding ions. 

Our main finding is that, as expected, tensile strain 
lowers the oxygen vacancy formation energy in CaMnOs. 
This shows that the material will respond to tensile strain 
not only by by increasing its bond lengths and modify- 
ing its octahedral rotation angles, but also by reducing 
its oxygen stoichiometry by creation of oxygen vacan- 
cies. In addition, and perhaps unexpectedly, we find that 
strain modifies the oxygen vacancy formation energy of 
inequivalent oxygen sites markedly differently. This find- 
ing implies a possible strain control of oxygen vaca ncy 
ordering, which could enable new functionalities^ 2 ^. 



II. COMPUTATIONAL DETAILS 

Density functional theory (DFT) calculations were 
performed using the VASP code^MSJfD with the PBEsol 
functional 3 - 3 -' and wave functions expanded in plane waves 
up to a kinetic energ y cuto ff of 650 eV. The projector aug- 
mented wave 

(PAW^PMl method was used to describe 
electron-core interactions with Ca(3s, 3p, 4s), Mn(3p, 3d, 
4s) and 0(2s, 2p) shells treated as valence electrons. The 
Hubbard U correctional was applied to Mn 3d electrons, 
using a U value of 3.0 eV, which was found to repro- 
duce the experimentally determined density of states^ 
and was recently shown to yield the correct stability of 
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the different magnetic phases of CaMnOsP^ as well as the 
redox energetics of binary manganese oxides^. We used 
the ground-state G-type antiferromagnetic order in con- 
junction with a 2 x 2 x 2 supercell of the 5 atom perovskite 
unit-cell. Reciprocal space integration was performed on 
a 4 x 4 x 4 T-centered mesh. Structures were set up within 
the desired Pnma space-group and subsequently relaxed 
until forces on the ions converged to below eV/A. In 
our calculations we strained CaMnC>3 in the plane of the 
equal a and c axes, while relaxing the &-axis length and 
all internal coordinates. During calculations containing 
Vo the lattice parameters were kept fixed respectively at 
the strained a = c and relaxed b lattice parameter, while 
relaxing only internal coordinates. 

III. RESULTS & DISCUSSION 
A. Stoichiometric CaMnO;i 

We first look at strain induced changes in the stoichio- 
metric CaMn03 structure, characterised by the volume, 
lattice parameters, bond lengths and octahedral rotation 
angles shown in Fig. [3] We notice first in Fig. |3ja), 
that strain leads to a net volume change in CaMn03, 
meaning that the material does not have a perfect Pois- 
son ratio. This is also reflected in the lattice parameters 
shown in Fig. [3jb). While the a and c lattice parameters 
change linearly as a function of the imposed strain, the 
b lattice parameter expands more strongly in the com- 
pressive range than it shrinks in the tensile range. It's 
rate of change is however about a factor two smaller than 
required for constant volume, leading to the observed vol- 
ume decrease with increasing compression. 

As schematically shown in in Fig. [l] in the stoichio- 
metric structure, these changes in lattice parameter are 
accommodated by changes in bond-lengths and octahe- 
dral rotation angles. Indeed in Fig. |3jc) and (d) we see 
that both the bond lengths and tilt angles change as a 
function of bi-axial strain. The bond lengths are how- 
ever relatively insensitive to strain (an order of magni- 
tude smaller change than the lattice parameters, see Fig. 
|3jb)), meaning that for stoichiometric Pnma CaMn03 
bi-axial strain is primarily accommodated by altering the 
octahedral tilt angles, Fig. J3jd) . As expected, the octa- 
hedral tilt angles increase with compression of the respec- 
tive axis, and decrease when the axis becomes longer. 

B. Tensile strain 

In the orthorhombic CaMn03 structure there are two 
inequivalent oxygen positions as shown in Fig. [2] We 
label these IP and OP, corresponding to the broken Mn- 
O bonds lying in- or out-of the plane of epitaxial strain 
respectively, but note that the sites are already inequiva- 
lent in the unstrained Pnma structure. In Fig. [4]we show 
our calculated formation energies for these two oxygen 
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FIG. 3. Evolution of structural properties with biaxial strain 
in stoichiometric CaMnOs: a) relaxed 8- formula- unit super- 
cell volume (a material with an ideal Poisson ratio would have 
a constant volume with strain), b) lattice parameters (the 
dashed line shows the b-axis parameter for a material with 
ideal poisson ratio), c) Mn-0 bond lengths and d) octahedral 
rotation angles (defined as angles between cartesian axes and 
Mn-0 bonds). 



vacancies as a function of the applied strain at an oxy- 
gen chemical potential of -5 eV, corresponding to typical 
growth conditions under air. In the inset in Fig. [4] we 
show the Vo formation energy in the 2% tensile strained 
structure as a function of the oxygen chemical potential 
throughout the whole range of CaMn03 thermodynamic 
stability evaluated by computing the Ca-Mn-0 ternary 
phase diagram^. As is well known, the oxygen chem- 
ical potential strongly affects the absolute value of the 
formation energy and hence the equilibrium oxygen va- 
cancy concentration. 



1. Change in stoichiometry 

As expected from the known chemical expansion 
caused by oxygen vacancies, we find that tensile strain 
lowers the formation energy for oxygen vacancies in 
CaMn03. The formation energy decreases considerably 
by ^0.4eV corresponding to ~20% of the bulk value for 
formation of the IP vacancy at 4% strain. 

In Fig. [5ja) we show the calculated equilibrium change 
in stoichiometry as a function of temperature for tensile 
strain values of and 4% and oxygen partial pressures 
corresponding to growth in air (0.21) and under nitro- 
gen (10 -6 ). The absolute oxygen stoichiometry 3-<5 was 
estimated from reaction (n}, the equilibrium constant of 
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FIG. 4. Vo formation energies for the two inequivalent Vo 
positions in the Pnma tilted CaMnOa structure at an oxygen 
chemical potential of -5 eV, corresponding to typical growth 
conditions under air. The inset shows the formation energy 
of the IP Vo at 2% tensile strain as a function of the oxygen 
chemical potential. The OP Vo very closely tracks this curve 
and is omitted for clarity. 



which can be written 

[Mn»+] 2 [Vo]pg 
[Mn 4 +] 2 [0 2 ^] 



K 



exp 



AGf 



(2) 



Here Gf is the Gibbs energy of reaction 
reaction ([I]) and charge neutrality we get 



( TJ . From 
Vo] = 5, 

[Mn 3+ ] ^25, [Mn 4 +] =1-25 and [0 2 ~] = 3 - 5. 
Equation ^ was solved for 5 using AGf = AHf — TAS 
where AHf is the enthalpy of oxygen vacancy formation 
and AS the entropy of reaction Q, respectively, where 
the latter is 95 JK _1 mol _1 taken fromSI. 

At a typical growth temperature of 1000 K (indicated 
by the vertical dashed line in Fig. [5]), tensile strain of 4% 
reduces the equilibrium oxygen stoichiometry from 2.995 
to 2.978 in air and from 2.964 to 2.859 in N2 atmosphere 
(po 2 = 10~ 6 atm.), as illustrated in Fig. [5] The oxygen 
stoichiometry of the supercell used in our calculations is 
3-5 = 2.875 for comparison. A reduction of 0.4 eV of the 
Vo formation energy has an impact on 5 equivalent to 
raising the temperature 150-200 K. and will drastically 
alter the equilibrium oxygen content. 



2. Vacancy ordering 

While both vacancy types lower their formation en- 
ergy, we find that the energy lowering is stronger for the 
IP vacancy than for the OP, consistent with the fact that 
tensile strain expands the in-plane lattice parameter and 
leads to a slight contraction of the out-of plane dimen- 
sions (see Fig. [3] (b)), however still with a net increase 
in volume as shown in Fig. (BTa). This suggests that, in 
addition to tuning the oxygen vacancy formation energy, 
tensile bi-axial strain could be used to influence vacancy 
site preference and thus induce vacancy ordering. 
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FIG. 5. a) Modelled oxygen stoichiometry of CaMn03_f. The 
3-5 curves are calculated from equation see text for ex- 
planation, b) Temperature dependent Vo ordering expressed 
as the ratio of IP/OP (Fig. [5} for +2% and +4% tensile 
strain (Fig. j4j). The typical growth temperature of 1000K is 
indicated by the vertical dashed line. 



In Fig. (5^b) we show our calculated equilibrium ratio 
of IP to OP Vo as a function of temperature for tensile 
strain values of 2 and 4%. The ratio of IP to OP va- 
cancies under tensile strain was calculated from the shift 
in formation energies shown in Fig. [4] Under +4% ten- 
sile strain, for example, the IP Vo formation energy is 
0.36 eV lower than that of the OP Vo- At the typical 
growth temperature of 1000KP^ (indicated by the vertical 
dashed line in Fig. [5]) this energy difference corresponds 
to an equilibrium ratio of IP to OP Vo = ~73. How- 
ever, oxygen vacan cies can remain mobile to much lower 
temperatures^lMI. Equilibrium vacancy concentrations 
can thus be reached at significantly lower temperatures, 
especially in thin films where the diffusion distances are 
shortP. At 300 K the equilibrium ratio of IP to OP va- 
cancies is more than 10 6 , (see Fig. [5](b) demonstrating 
the possibility of engineering oxygen vacancy ordering 
by epitaxial strain. Oxygen vacancy ordering induced 
by epitaxial strain modification of the vacancy forma- 
tion energy can in principle be applied to any oxide with 
non-equivalent oxygen sites. This can pave the way to- 
wards new interface functionalities arising from orbital 
ordering (Fig. [6] (d)), localized electric dipoles/fields^ 
or enhanced interface electronic, ionic or mixed conduc- 
tivity. 



C. Compressive strain 

In spite of the fact that the volume of the cell de- 
creases monotonically from tensile to compressive strain 



5 



as shown in Fig. [3^a), we find a negligible change in 
the energy of formation of either vacancy on compressive 
strain. While our findings for tensile strain are consis- 
tent with the simple arguments based on changes in Mn 
ion valence and volume presented in the introduction, 
the compressive-strain response is fundamentally differ- 
ent. Next we analyze the changes that strain induces in 
the electronic structure to understand the origin of this 
behavior. 

In Fig. [61(a) we show our calculated electronic density 
of states (DOS) for unstrained CaMn03 with one out of 
24 O atoms missing. We see that IP Vo formation pro- 
duces an extra peak within the band-gap of the density of 
states. The OP Vo produces a peak at nearly the same 
position and has hence been omitted for clarity. The 
IP Vo peak accommodates the two charge compensating 
electrons, which are primarily located in the formerly un- 
occupied z 2 -type orbitals of the two Mn atoms adjacent 
to the oxygen vacancy, as shown in the calculated charge 
density in Fig. |6](d). Localisation of these electrons is 
a result of the removal of the anti-bonding a* overlap 
between the Mn e g and the O p orbitals along the axis 
of the broken bond, which lowers the energies of these e g 
orbitals. 

We also observe that both vacancies will cause an up- 
ward shift of about ~0.1 eV of the valence band edge, 
which consists of Mn ti g orbitals. For the IP structure 
the increase stems from the xy and xz orbitals, whereas 
in the OP structure the xy and yz orbitals move up in 
energy. This reflects the loss in Mn-0 7r bonding overlap 
with the bond broken by IP Vo formation lying along 
the a (x) axis, whereas that broken by OP Vo lies along 
the b (y) axis. 

In Fig. [6] (b) and (c) we look at the changes in band 
edge and gap state in the DOS of the IP and OP structure 
compared to the stoichiometric case. We can see that the 
valence band edge for the IP Vo is lowered by ~ 0.05 eV, 
whereas for the OP Vo is raised by an equal amount. We 
can again rationalize these effects in terms of structural 
changes. In the stoichiometric and IP case, the valence 
band-edge consists primarily of xz orbitals lying in the 
compressed plane. Creation of an Vo in this plane can 
to some extent counteract the increase in anti-bonding 
overlap induced by compression and hence lead to a slight 
net lowering of the band edge energy in the IP case. This 
is not the case for the OP Vo, for which moreover the 
xy and yz orbitals increase in energy due to loss of Mn-0 
bonding overlap, resulting in the observed net upwards 
shift. 

The opposite trend is observed for the gap state, which 
is lowered in energy by ~ 0.1 eV for the OP case, whereas 
in the IP case we observe a broadening with the lower 
edge dropping by ~ 0.05 eV and the upper edge remaining 
more or less constant. The downwards shift of the OP 
gap state is due to decreased a* overlap since the broken 
Mn-0 bonds lie along the stretched b axis. Conversely 
the broadening for the IP case stems from an increase of 
this overlap due to compression of the Mn-0 bonds lying 
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FIG. 6. Mn d partial density of states shown for (a) the 
unstrained IP Vo case, (b) the IP Vo at 4% compressive 
strain and (c) the OP Vo at 4% compressive strain. The 
DOS of the stoichiometric material is shown in grey for the 
corresponding strain. The energy scales of the stoichiometric 
DOS have been aligned with respect to the defective DOS at 
a Mn p core-state about 50 eV below the Fermi energy. The 
origin in energy is set at the Fermi energy of the defective 
cell. Panel (d) shows the charge density associated with the 
occupied e g states upon IP Vo formation. The dashed red 
lines and circle indicate the bonds and position of the missing 
O atom. Yellow and light blue isosurfaces at 8 • 10 -3 eA -3 
show the up and down spin density respectively. Color code: 
Ca: blue, Mn: purple, O: red. 



in the compressively strained plane. 

In the compressive regime we thus have counterbalanc- 
ing effects, leading to very similar formation energies for 
the two vacancies, which are moreover almost indepen- 
dent of strain. 

The above results indicate that crystal field effects 
dominate changes in Vo formation energies under com- 
pression. This can be further assessed by a decomposition 
of the Vo formation energy into its components. In Fig. 
[7] we show that in the tensile regime the band energy re- 
lated part (Fig. |7|i) alone would predict the OP Vo to 
be lower in energy. In order to recover the observation 
in Fig. [4] that the IP Vo is more stable, the remaining 
- mainly electrostatic - interactions (Fig. [7|d) have to 
dominate in the tensile range. The marked difference be- 
tween IP and OP under tension is a result of the drastic 
reduction in repulsive electron-electron interaction due 
to stretching of the broken bond axis. Conversely in the 
compressive range in Fig. [7] we can see that the observed 
behaviour of the total formation energy (Fig. [4]) stems 
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FIG. 7. The total formation energy shown in Fig. [^decom- 
posed into contributions from a) the band energy and b) the 
remaining interactions. 



mainly from changes in the band energy due to crystal 
field effects. 

In our calculations relaxation of the internal coordi- 
nates allowed both octahedral rotation angles and bond 
lengths to vary as a function of strain. As mentioned 
earlier (Fig. [3]). changes in octahedral rotation angles 
account for a large part of strain accommodation in the 
Pnma structure. Therefore crystal field-effects due to 
bond-length changes should be rather weak. In order 
to asses the limit of strong crystal-field effects, when no 
changes in tilt pattern exist , we have performed the same 
analysis on the ideal cubic structure. Under tensile strain 
both vacancies show a similar although less marked be- 
havior than in the Pnma structure. In the compressive 
range however we find the OP Vo formation energy to 
continuously increase whereas the IP Vo shows a marked 
decrease. This is the result of a pronounced broadening 
of the e g gap states in the cubic structure as the strain is 
accommodated by changes in bond-length (overlap) only. 
These observations confirm again that under compression 
crystal-field effects dominate whereas electrostatics are at 
the origin of the observed behaviour under tensile strain. 

IV. SUMMARY & CONCLUSIONS 

In the present work we have investigated how the oxy- 
gen vacancy formation energy is modified by bi-axial 
strain from 4% compressive to 4% tensile strain. Our 
DFT calculations show that oxygen vacancy energetics 
are strongly modified by strain, particularly by tensile 
strain which was found to favor the formation of oxy- 
gen vacancies. Compared to the unstrained material, the 
formation energy of oxygen vacancies is reduced by up 
to 0.4 eV under 4% tensile strain, significantly reduc- 



ing the equilibrium oxygen content in films under tensile 
strain. The origin of this reduction is mainly electro- 
static, as tensile strain reduces the electron-electron re- 
pulsion along the broken Mn-O-Mn bond axis. It thus fol- 
lows that CaMnOa can accommodate strain b y increasing 
its concentration of oxygen vacancie d 10 * 21 * 44 ! 

Compressive strain on the other hand was found to af- 
fect Mn d levels in a more subtle way with nearly no 
net effect on oxygen vacancy formation energies. Al- 
though not explicitly treated in this work, we expect the 
cation vacancy concentration in oxygen hyperstoichio- 
metric perovskites, like LaMn03 +( 5, to be correspond- 
ingly sensitive to compressive epitaxial strain^. In such 
materials the hyperstoichiomctry is a result not of oxygen 
excess but cation deficiency. 

This stabilization of oxygen vacancies by tensile strain 
is the reverse effect of the well established chemical ex- 
pansion phenomenon where oxygen vacancies induce a 
unit cell expansion^! Our findings have profound impli- 
cations for the properties, growth and thermal annealing 
of oxide thin films under tensile strain. Epitaxial strain, 
like temperature and partial pressure of oxygen (po, ) , is 
a thermodynamic variable affecting the equilibrium oxy- 
gen stoichiometry. 

The formation of oxygen vacancies must be charge 
compensated, either by reducing the formal oxidation 
state of cations (e.g. Mn 4+ to Mn 3+ in the present work) 
or more unfavorably by donating electrons into the high 
energy conduction band (e.g. LaAlOjpQl). The formation 
of oxygen vacancies is therefore likely to be enhanced by 
tensile strain in any oxide with transition metal cations 
with multiple oxidation states such as Mn, Fe, Co, Ni 
or even Cu. The coupling between strain and oxygen 
vacancy formation energy is moreover expected to be 
most pronounced in materials containing transition metal 
cations with localized d-electroniPSEZI 

Our calculations also show that tensile strain affects 
the formation energy of the two kinds of non-equivalent 
oxygen vacancies in CaMnOa with Pnma perovskite 
structure in a markedly different way. Under 4% ten- 
sile strain the formation energy of oxygen vacancies on 
the in-plane site is 0.36 eV lower than on the out-of-plane 
site. Depending on the temperature this large energy dif- 
ference may result in orders of magnitude differences in 
vacancy concentrations on the two sites and paves the 
way for tensilc-strain-induced oxygen vacancy ordering. 
This possibility of engineering epitaxial heterostructures 
with oxygen vacancy ordering should apply to any oxide 
with non-equivalent oxygen sites. 
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